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Abstract: A coupling kinetics model is developed to simulate the release and transport of 
landfill leachate pollutants in a deformable municipal solid waste landfill by taking into 
account of landfill settlement, seepage of leachate water, hydrolyse of insoluble and 
degradable organic pollutants in solid phase, biodegradation of soluble and degradable 
organic pollutants in solid phase and aqueous one, growth of aerobic and anaerobic 
microorganism, and consumption of dissolved oxygen. The release and transport of organic 
pollutants and microorganisms in landfills in the process of landfill settlement was 
simulated by considering no hydraulic effect. Simulation results demonstrated that the 
interaction between landfill settlement and the release, transport and biodegradation of 
landfill leachate pollutants was significant. Porosity and saturated hydraulic conductivity 
were not constants because of the landfill settlement, which affected the release, transport 
and biodegradation of landfill leachate pollutants, and furthermore acted on the landfill 
settlement. The simulation results accorded with the practical situation, which preliminarily 
verified the reliability of the mathematical model and the numerical program in this paper. 

Keywords: municipal solid waste; landfill leachate pollutants; aerobic and anaerobic 
degradation; dissolved oxygen; settlement; coupling kinetics model 
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1. Introduction 

The release and transport of landfill leachate is a complex process, affected by landfill settlement, 
fluid movement, biodegradation and temperature changes, so a complete model which describes the 
release and transport of landfill leachate pollutants must contain mechanical, hydraulic, gas transport, 
temperature and biodegradation models, but it is almost impossible to realize a five-field coupling 
simulation, so this is often simplified to two or three field coupling model. 

Many researchers have studied the multi-field coupling problems of landfill leachate transport, and 
new models have been developed based on more detailed mathematical descriptions of the landfill and 
incorporating other aspects of interest apart from hydrology, such as the biological and physical- 
chemical degradation and settlement. Demirekler et al. [1] developed a three-dimensional 
mathematical model to estimate the quality and quantity of the landfill leachate produced. The effect of 
overburden stress was considered. Lobo et al. [2,3] has reported the first version of MODUELO. The 
development of the second version and the Meruelo Landfill (Spain) simulation results were presented 
in Lobo et al. [4-6]. Chanthikul et al. [7] developed a mathematical model of BOD5 concentration 
without and with leachate recirculation. Durmusoglu et al. [8] developed a one-dimensional multiphase 
numerical model to simulate the vertical settlement involving liquid and gas flows in a deformable MSW 
landfill. McDougall [9] developed a hydro-bio-mechanical model for settlement and other behavior in 
land filled waste. Fellner and Brunner [10] established a 2-dimensional 2-domain model for simulating 
the leachate generation from MSW landfills. A flow field consisting of a vertical path (channel domain) 
surrounded by the waste mass is defined using the software HYDRUS-2D. One-dimensional advection- 
dispersion transport modeling was conducted as a conceptual approach for the estimation of the transport 
parameters of fourteen different phenolic compounds and three different inorganic contaminants 
migrating downward through the several liner systems in Gamze et al. [11]. Gaetano et al. [12] presents 
a ID mathematical model for the simulation of the percolation fluxes throughout a landfill for MSW, 
which considered the landfill divided in several layers evaluating the inflow to and outflow from each 
layer as well as the continuous moisture distribution. But, there are still one or more defects in most 
models as follows: (a) The pollutants were considered as a single solute, whether it is soluble or 
insoluble, degradable or non-degradable was not definitely differentiated; (b) Solid-phase pollutants 
and its dissolution to aqueous phase weren't taken into account; (c) The effect of oxygen on 
biodegradation was neglected and the transition from aerobic degradation to anaerobic one wasn't 
taken into account; (d) the effect of settlement on porosity, saturated hydraulic conductivity and the 
transport of pollutant was neglected. 

A coupling kinetic model was developed to simulate the release and transport of leachate pollutants 
in a deformable MSW landfill taking into account of hydrolyse and dissolution of solid-phase 
pollutants, oxygen consumption and transition of aqueous-phase pollutant biodegradation from 
anaerobic stage to aerobic one, and other behaviors such as convection and hydrodynamic dispersion, 
adsorption/desorption and growth of microorganism. A case study was given by considering none 
hydraulic action for studying the change law of water quality and quantity, which preliminarily 
verified the reliability of the mathematical model by comprising with the practical situation. 
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2. Mathematical Model 

2.1. Basic Assumptions 

The release and transport of organic pollutants in landfill is a complicated process which is 
accompanied by physical behavior and chemical and microbial reactions. It can be barely described by 
a completely correct model. The development of the simulation model must be based on some suitable 
assumptions. The assumptions of the models in this study are as follows: (a) Landfill gas is released 
rapidly after generation, so the landfill leachate transport is considered as a single phase flow; 
(b) MSW particles are incompressible, but degradable; (c) The simulated landfill was taken as a 
biochemical reactor. Organics transport and transform under a series of physical, chemical and 
biological actions, such as convection and hydrodynamic dispersion, hydrolyse, dissolution, 
adsorption/desorption and biodegradation; (d) Density and viscosity coefficient of landfill lecheate are 
constants. 

2.2. Landfill Settlement Model 

2.2.1. Mass-Conservation Equation 

Based on the mass conservation principle the mass-conservation equation of solid phase is: 

d(l-g>)jp, dQ.-q>)p s v; 

dt dx t Us CI) 

where x, is the space coordinates [L]; r is the current time [T]; p s is the solid phase density [ML ]; 
<p is the porosity; v. is the velocity of solid phase [LT 1 ]; Q s is the source/sink term, which is caused 

—3 — 1 

by the degradation of solid waste, the release of inner source water etc. [ML T ]. 

2.2.2. Mechanical Model 

The Merchant model was used to simulate landfill settlement. It was constructed by a Hooke 
elastomer and a Kelvin model in series. Kelvin model was constructed by a Hooke elastomer and a 
Newton viscosity mode in parallel. The creep equation is: 



where e i} is the deviator strain tensor [LL Sy is the deviator stress tensor [ML l T 2 ]; K ] is the 
viscosity coefficient [ML T ]; G 0 and G l are the shear modulus for Hooke elastomer and Kelvin 
model, respectively [ML T ]. 




(2) 



In addition: 



e =£..—£ 8 



(3) 



s.. = a '..-a' 8.. 



(4) 
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So the effective stress can be written as: 



(5) 



In Equations (3)-(5), s tj is the strain tensor [LL ]; e n = {e x +e y +s z )l3 is the mean strain [LL ]; cr^ is 
the effective stress tensor [ML T ]; a' m =(a' x + a' y + a' z )/3 is the mean effective stress [ML T ]; A" is 

-W2- 



the bulk modulus [ML T ]; S v is the Kronecher symbol; and: 



2G 0 2G 1 



l-e 



(6) 



Furthermore, the effective stress principle can be described by: 



a = <t' + S v8 = — s + 



where a tj is the total stress [ML ! T 2 ]; S w is the liquid saturation [LL 3 ]; p is the liquid pressure 
[ML _1 T~ 2 ]. 

Geometric equation and stress equilibrium equation are: 

1 



(3K--)e m + S w p 
A 



(7) 



•./;-(> 



(8) 



(9) 



The stress equilibrium equation represented with displacement can be obtained by plugging 
Equation (7) and Equation (8) to Equation (9): 



(3^ + ^+^+^ + /,=0 



The velocity of solid phase is: 



(10) 



(11) 



Equation (1) and Equation (7) (or Equation (10)), Equation (8) and Equation (11) are the basic 
equations of landfill settlement model. The liquid pressure p was contained in it, so the hydraulic 
model must be developed for obtaining p . 

2.3. Hydraulic Model 



Based on the mass conservation principle the continuity equation of aqueous phase is: 



dt 



dx. 



Qw 



(12) 



where p w is the liquid phase density [ML 3 ]; g H is the source/sink term [ML 3 T ']; v\\ =vj+v^, is the 
absolute velocity of aqueous phase [13] [LT 1 ]; and vj is the relative velocity of aqueous phase to the 
solid phase [LT 1 ]. 



Int. J. Environ. Res. Public Health 2012, 9 



3441 



During settlement, the solid particles as well as the liquid move simultaneously. Hence, it is 
necessary to state Darcy's law relative to solids movement. That is: 



6vT = 



dh 

kijk r (h) \-k iz k r (h) 

dxj 



(13) 



where 6(h) = 0S W is the volumetric moisture content [L 3 L 3 ]; his the hydraulic pressure head [L]; k tj is 
the saturated hydraulic conductivity tensor [LT 1 ]; k r (h) is the relative permeability. 
VG function [14] is used to describe the water retention curve: 



0(h) = 0 r + (0 s -0 r )(l + \ah\ n f m (14) 



where a ,n ,m = \-\ln are parameters; 6 r and 6 S are the residual and saturated volumetric moisture 
contents, respectively. 

So the relative permeability can be written as: 

k r (h) = Sl l2 [l-(l-Sl"Tf (15) 

where S e = ( 6 - 6 r ) I ( 9 S - 9 r ) is the effective saturation. 
ky and Q w [15] are calculated by: 

i 

k..=A- (e 8 ^ - e yB ) (i = j) or h. = 0 (/ * j) ( 1 6) 

Q K = ( Q 0 + A, • e~' IB > + A, • e-" B > ) • (1 - <p)p s ( \ 7) 

where A , B ,Q Q , A l , A 2 , 6j and fi 2 are parameters; p s is the particle density of MSW [ML ]. 

2.4. Pollutant Release and Transport Model 
2.4.1. Conceptual Framework 

Organic biodegradation in landfills can be divided into two stages: (1) aerobic biodegradation and 
(2) anaerobic biodegradation. The first one always occurs in the initial landfill stage, and it can be also 
divided into two stages: (1) the hydrolysis stage of insoluble macromolecular organics to soluble and 
small molecular ones and (2) the biodegradation of soluble organics to H2O and CO2, etc. When the 
oxygen is consumed, biodegradation enters the anaerobic stage. In this stage, macromolecular organics 
are hydrolyzed to small molecular ones, and then decomposed to CH 4 and H 2 0 by anaerobic 
microorganisms after the acidification process. 

Based on above biodegradation process, organic pollutants in landfill can be classified as insoluble 
and degradable ones (IDS), soluble and degradable ones (SDS), and adsorbed ones (AS) in solid phase; 
and soluble and degradable ones in aqueous phase (SDA). Microorganism includes aerobic and 
anaerobic ones in aqueous phase (AM and ANM) and hydrolysis ones in solid phase (MS). The model 
which describes the pollutant release and transport in landfill can be developed by using the mass 
conservation principle, including hydrolysis of IDS, dissolution and biodegradation of SDS, 
adsorption/desorption and aerobic and anaerobic biodegradation of SDA; growth and death of AM, 
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ANM and MS, and consumption of dissolved oxygen (DO). The biodegradation process of organics 
was shown in Figure 1 . 

Figure 1. The biodegradation process of organics in landfills. 
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if Anaerobic biodegradation 
I "3 
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Aerobic biodegradation 



Methane, etc. final product 



H 2 0, C0 2 , etc. final product 



2.4.2. Hydrolysis of IDS 

The hydrolysis of insoluble macromolecular organics (IDS) to soluble and small molecular ones can 
be described as first order reaction: 



nsi is ni 

1<H - -&SiH>J 



(18) 



where R% is the hydrolysis rate of IDS [MM T ]; 5' is the concentration of IDS [MM ] and K sm 
is the hydrolysis constant [T ]. 

2.4.3. Dissolution of SDS 

The dissolution of SDS is closely related to the water content and the pollutant concentrations in 
solid and aqueous phase. It is described by [16]: 



rt d 

>Sd ^ - 



Rs - ~{~Zj) m KsdS (Cmax - C" )# 



(19) 



where Rs" is the dissolution rate of SDS [MM ! T S" and So are the concentrations of SDS at time t 
and initial time, respectively [MM -1 ]; C c is the concentration of SDA [ML -3 ]; Cmax is the maximum 

—3 —3 — 1 

concentration of SDA [ML ];K S dsis the dissolution rate constant [ML T ]; m is the dissolution 
coefficient. 

2.4.4. Biodegradation 

The decomposition and stabilization of MSW in landfill is essentially a microbial metabolic process. 
The depletion of the substrate and microorganism growth can be described by Monod kinetics [17], 
hence for MS accumulation: 



Rd=^ 



max SdD 



KsdD + S Cl 



(20) 
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where R s D m is the growth rate of MS [MM -1 T -1 ]; S m is the concentration of MS [MM -1 ]; /J maxS dD is the 
maximum specific growth rates for MS [T 1 ]; K sdD is the half saturation constant for SDS [MM -1 ]. 

The depletion rate of the substrate is directly related to MS accumulation through a cell/substrate 
yield coefficient Y Sm/sd : 



11 S d 

Rd = ~- Rd" 1 = ~- AnaxSdD " — —7 ■ S" (21) 



where Rd is the depletion rate of SDS [MM T ]; the Ysm/sd is the stoichiometric yield coefficient for 
MS (biomass produced per unit amount of electron donor utilized) [MM - ]. 

When the dissolved oxygen (DO) exists, and its concentration is low, the cell growth rate for AM 
and ANM can be represented by the following double Monod models [18]: 



CC s~lO 
,,„ _ f , ^ 

"max CcO v f~ 
K CcO + C 

f r^o \ 



» ~MmsxCcO' v * C c Y J. C 



jyCma _ C_ 

D A'maxCeA 



K CcA + C 



C° 

1 



K CoA + C " J 



■ c 



ma 



(23) 



where i?£ mo and R 0 ™ are the cell growth rates for AM and ANM, respectively [ML -3 T -1 ]; C mo and 

C™ are the concentrations of aerobic microorganism and anaerobic microorganism, respectively 

[ML -3 ]; C° is the DO concentration [ML -3 ]; // maxCe0 and // maxCcA are the maximum specific growth 

rates for aerobic and anaerobic microorganism, respectively [T -1 ]; K Cc0 and K CcA are the half 

saturation constant for aerobic and anaerobic microorganism, respectively [ML ]; K Co0 andK CoA are 

the half saturation constants for DO [ML ]. 

The depletion rates of the substrate and DO in aqueous phase are directly related to AM and ANM 
accumulation through cell/substrate yield coefficients Y Cma/Cc ,Y Cma/Cc and Y ColCc : 

1 1 C c C° 

TfCc _ nCmo 1 ^ f«io /n a\ 

y ' D y '/VaxCcO'^ , C c ' K + ' ' 

1 CmolCc 1 CmolCc ^CcO + ^ A CoO + ^ 



E>Ce ; r> ; 

A T , ' ^CmaD' ' A*ir 



A CmaZ) ^ "max CcA j_ /"~" c 

*CmalCc *CmalCc ^CcA +( ^ 



1 

V K CoA +C° J 



■C" 



(25) 



1 c c c° 

pCo _ y nCc _ _y * .. //-«/'\ 

^CmolCc K CcO + ^ A CoO + ^ 

where R^ and i?^ c are the aerobic and anaerobic degradation rates of SDA, respectively [ML -3 T -1 ]; 

Rq" is the consumption rate for DO [ML -3 T -1 ]; Y Cmo/Cc and Y CmalCc are the stoichiometric yield 

coefficients for aerobic microorganism and anaerobic microorganism, respectively (biomass produced 
per unit amount of electron donor utilized); Y Co/Cc is the consumption coefficient of DO (oxygen 

consumed per unit amount of SDA). 

The MS, AM and ANM decay are given by: 

R s E m =-K SmE S m (27) 
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R c E mo =-K CmnE C mo (28) 



R c E ma =-K CmaE C' na (29) 



where R s E m , R E mo and Bg" are the endogenous cell death or decay rates of MS, AM and ANM, 



—3 — 1 

respectively [ML T ]; K SmE , K CmoE and K CmaE arQ the endogenous cell death or decay coefficients of 
MS, AM and ANM, respectively [T 1 ]. 

2.4.5. Adsorption/Desorption of SDA 

Langmuir adsorption model can describe the adsorption behavior of the pollutants in MSW well [19]. 
So the adsorption rate is described by: 



RZ=a<KS s m7 £—-n (30) 
L +K d 

where Rp is the adsorption rate [MM -1 ]; S s and 5* are the adsorption concentration and maximum 
adsorption concentration of IDS [MM \,K d is the equilibrium sorption constant [ML ]; a is the first 
order adsorption/desorption rate constant [T -1 ]. 

2.4.6. Governing Equations 

Based on the mass conservation principle and considering landfill settlement and hydrolysis of 
macromolecular organics, the governing equation for IDS can be described by: 

dp,(i-<p)s' dv;p,(i-<p)s l 



dt Bxj 



-p s (\-(p)K SiH S' (31) 



The SDS governing equation considering landfill settlement, hydrolysis of IDS and dissolution and 
biodegradation can be given by: 

dp s (l-(p)S d dv!p,(l-<p)S d , n \( S \™ v , r c r c^ a 

= + p s (1 - <p)Ys dH K SiH S - p s (1 - q>)(—j) K sds (C max - C )0 

dt dXi S 0 

1 s d ( 32 ) 

"A (1 " <P) ~ M™*sdD ■ — -j S 



The MS governing equation considering landfill settlement, growth and decay of MS is described as: 

d Ps (i-9)s m _ dt P ,<x-<p)sr s d 

+Ps(L-<P)-MmmSdD-— ~J ^ ~ AU~ <P) K SmE?> (33) 

dt 8^ K sdD +S v ' 

The AS governing equation is: 

+ PsKi-(P)atA?m— c — - — J ) (34) 

dt dXi C^+K d 

The SDA transport model considering convection and hydrodynamic dispersion, dissolution of 
SDS, aerobic and anaerobic degradation and adsorption/desorption is given by: 
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dec _ a ac c 

Sf cbc. '' dx. 



- + Op, (1 - ^ (C m „ - C c ) - 0- — — 



c° 



• c 



a 1 C 

_t/ ' /"max 



K CoA + C° 



(35) 



C™- A (l-p)a0(S* 



-5 s ) 



The equations for AM and ANM by considering convection and hydrodynamic dispersion, growth 
and decay are given by: 



d6C mo _ d 
dt dxj 

dec"" e 



dC 

ox 



d6v7C mo C 

-+ C//J meai CcO 



C° 



i J 



dt dx/ 



9D tj 



8C m 



dx 



j J 



dxi 



d0v?C' r - 
dx; 



KccO + C Kr„n + C 



Cmo pjy- 



CcA ' 



KcaA + C C 



CoO 



c 



1— 

V Kcoa + C 



OJ^CmaE^ 



(36) 
(37) 



The governing equation for DO is: 



d9C" d 



dt dx. 



0D ; 



dC" 



" dx 



J J 



dOv*;c" 

dx. 



-er, 



ColCc 



Mmax Co 



L Cmo/Cc 



K CcO +CC K CoO +C " 



■C" 



(38) 



2.4.7. Numerical Solution Method 



The Merchant model was obtained by Lagrangian description. Hydraulic model and pollutant 
release and transport model were obtained by an Eulerian description. Total settlement in untreated 
landfilled MSW has been estimated to range between 25% and 50% of initial fill height [20]. So the 
upper boundaries of fluid and pollutant transport regions are obviously moving, and the small 
deformation assumption isn't suitable. The coupling of these two types of model may lead to the 
moving boundaries of Eulerian describing models, so the Arbitrary Lagrangian-Eulerian (ALE) 
method was used to the model solution for solving the moving boundary problem. Due to space 
limitations, the solution process can be seen in the authors' another work [15]. 

3. Results and Discussion 

An ideal landfill should have an effective seepage control system. After closure, it is in a relative 
independent state and can't be affected by the external hydraulic environment. In this study, the change 
law of the main physical and chemical variables and its effect on the pollutant transport was analyzed 
in an ideal landfill. The simulated landfill had a rectangular vertical section of 15 m in height and 20 m 
in width. All boundaries were impervious. The upper boundary can move freely, and the others are all 
fixed. The parameters are given in Table 1 . Physical and mechanical parameters were determined by 
testing, and biological parameters were obtained by parameter inversion. The results are shown in 
Sections 3.1-3.10. In Figures 1 and 2, z is the space coordinate of a certain particle at the initial time, 
that is, the corresponding space coordinate z at initial time was used to represent a certain particle. In 
the other figures, z is the space coordinates of a certain particle when the MSW was filled for 30 years, 
that is, the corresponding space coordinate z at 30 years after MSW was filled was used to represent a 
certain particle. 
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Table 1. Model parameters. 



Parameter 


Values 


Parameters 


Values 


Parameters 


Values 


Go 


1,320 kPa 


K CcO 


5.0 kg-m" 3 


/'max StID 


0.02 cT 1 


G, 


86.2 kPa 


K CM 


0.01 kg-m" 3 




0.0005 d" 1 


K s 


2.0xl0 5 cT 1 


K CcA 


1.2 kg-nT 3 




0.001 d" 1 


«o 


0m 


K CoA 


0.03 kg-m" 3 


c; n „ 


35,000 kg-nT 3 




0.52 


Ksi,,E 


0.0001 d" 1 


So 


0 kg-kg" 1 


e r 


0.21 


K-CmoE 


0.01 d" 1 


<; d 


0 kg-kg" 1 


a 


1.74 m" 1 


K-CmaE 


0.0002 d" 1 


So" 


0 kg-kg" 1 


it 


1.38 


YsdR 


0.5 


So 


0 kg-kg" 1 


P, 


864 kg-nT 3 




0.05 


c„ c 


0 kg-m" 3 


Ksm 


0.0006 d" 1 


V 

1 CmolCc 


9xl0" 5 


cr 


1.2xl0" 4 kg-nT 3 


Ksiis 


5xl0" 9 kg-kg" 


1 Y 


6x1 0~ 5 


C7 


2.5 xlO" 6 kg-m" 3 




5.0 kg-kg" 1 


V 

1 ColCc 


100 


r" 


8xl0" 3 kg-nT 3 




5.0 kg-kg" 1 


V 

1 Co/a- 









3.1. Displacement 

Figure 2 shows that the settlement occurred in almost 2 years. It's about 85% of total settlement. 
The total settlement was about 2.6 m, which was about 17.3% of initial fill height. The simulation 
results fitted well with the observed data from a similar landfill cell in Wuhan Jinkou landfill in China 
which was observed from 2001 to 2010, and it accorded with the reported settlement law [21,22]. 



Figure 2. Displacement change with time. 
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3.2. Porosity 

Figure 3 shows that porosity decreased at first and then increased with time due to the landfill 
settlement and organic biodegradation. The 85% settlement that occurred in 2 years led to the MSW 
compression and the decrease of porosity. After 2 years, the effect of biodegradation on porosity was 
more and more obvious. The organic biodegradation led to the reduction of solid mass, and the 
porosity presented an increasing trend. When the MSW was filled for 30 years, the porosities at top 
and bottom increased to 0.55 and 0.458 again, respectively. Meanwhile, the decrease of the MSW 
porosity at bottom due to the landfill settlement was more obvious; however the increase of the one at 
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top due to organic biodegradation was more obvious. In addition, the porosity presented an increasing 
trend with the fill height increasing. 

Figure 3. Porosity change with time. 
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3.3. Saturated Hydraulic Conductivity 

Figure 4 shows that the change of saturated hydraulic conductivity was similar to that of porosity, 
and the effect of displacement was significant. It decreased at first and then increased taking two years 
as a turning point. Increase was the main trend of Ks of the upper MSW, and the lower ones had a 
decreasing trend. The Ks at z = 0.73 m decreased from 0.8 m-day -1 , which was the initial value, to 
0.32 m-day -1 , and then tended to be stable; and the one at z = 12.4 m showed a small decrease at first, 
and then increased until the maximum value. 




3.4. Pressure Head 



It's seen from Figure 5 the effect of settlement on pressure head wasn't significant, although the 
water transport was closely related to landfill settlement, porosity and saturated hydraulic conductivity. 
The water mainly moved from top to bottom under gravity action and the upper pressure head 
decreased with time and the lower one increased by taking 5 m-6 m as separatrix. 
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The pressure head at the bottom reached 0.6 m when MSW was filled for 5 years, and increased 
gradually with time. When MSW was filled for 30 years, the MSW below 2.5 m was saturated, which 
was equivalent to 16.7% of the landfill height. Thus, although without the effect of groundwater and 
surface water invasion, the landfill leachate generated by MSW itself was large, and can't be neglected 
when designing the seepage control system and leachate treatment system. 

Figure 5. Pressure head change with time. 
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3.5. IDS 

Figure 6 shows that because it's assumed that the hydrolysis rate of IDS was only related to its own 
concentration, and the whole landfill cell has the same initial IDS concentration values, the calculated 
IDS concentration was only the decaying exponential function of time and independent of fill height 
and other parameters. 

Figure 6. IDS Concentration change with time. 
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3.6. SDS 



SDS concentration in this landfill increased at first and then decreased with time as seen in Figure 7. 
This change process was closely related to the dissolution and biodegradation of SDS and the 
hydrolysis of IDS. Because of the fast hydrolysis of IDS in the early period of the landfill, the SDS 
concentration increased rapidly. Over time, IDS concentration gradually decreased, and SDS 
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continuously dissolved from solid phase and was biodegraded, so the SDS concentration decreased 
year by year. The inflexion point was about the eleventh year. 

The difference of SDS concentrations between every height was small and is more obvious before 
the MSW was filled for 15 years. Because the pore water of lower waste was tending to be saturated 
gradually, and the dissolution of soluble organic pollutants was accelerated, there was a little 
difference between upper waste and lower one after the MSW was filled for 15 years. 

Figure 7. SDS concentration change with time. 
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3.7. MS 



The MS concentration increased year by year, and the difference between every height was small as 
seen in Figure 8. 

Figure 8. MS concentration change with time. 
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3.8. SDA 



Figure 9 shows the simulation results of SDA fitted well with the observed data from the Wuhan 
Jinkou landfill in China, which verified the accuracy of the model and parameters. These data were 
observed from 2001 to 2010. The suspended particles in water samples were filtered out and the SDA 
content was represented by chemical oxygen demand (COD) of the sample. SDA concentration 
increased at first and then decreased. At the early period, because of the higher SDS concentration and 
its dissolution to aqueous phase, the SDA concentration increased rapidly. Meanwhile, the water 
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content of the lower waste was relatively higher, which accelerated the dissolution of SDS, so the SDA 
concentration in lower waste was higher than the upper one. The rapid dissolution of SDS in the early 
period also led to the higher SDA concentration in upper waste than the lower one in the later period of 
landfill. 




It is seen from Figure 10 that the ANM concentration change was similar to that of MS in Figure 8, 
but it had a significant difference along height. The ANM concentration at z = 0.73 m was almost three 
times of the one at z = 12.40 m at utmost. This is mainly because the higher organic concentration in 
the lower waste provided sufficient nutrients for microorganisms and promoted their growth, while at 
later periods the SDA concentration in lower waste became lower, the ANM concentration presented a 
decreasing trend combining with its own decay. 

Figure 10. ANM concentration change with time. 
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3.10. AM and DO 

Figures 11 and 12 show that the change of DO was consistent with AM, which concentrations 
presented a rapid decreasing trend with time, and reached 0 when MSW filled for 200 days. The values 
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for the 1,001th day to 10,000th day were omitted in Figures 11 and 12 because they were the same as 
on the 1,000th day. 

Figure 11. AM concentration change with time. 
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Figure 12. DO concentration change with time. 
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4. Conclusions 



A coupling kinetic model of landfill leachate pollutant release and transport in the process of 
landfill settlement was developed, which contained three sub models-landfill settlement model, 
hydraulic model and pollutant release and transport model. Landfill settlement, convection and 
hydrodynamic dispersion of leachate, hydrolysis, dissolution, adsorption/desorption, biodegradation of 
pollutant and other behaviors were considered. The release and transport of pollutants and 
microorganisme in a landfill was simulated by considering no hydraulic action. The total settlement in 
this landfill cell was about 2.6 m, which was about 17.3% of initial height, and 85% almost occurred 
within 2 years. The simulation results fitted well with the observed data, and accorded with the 
reported settlement law. The changes of porosity and saturated hydraulic conductivity were closely 
related to settlement and biodegradation. They all presented a decreasing trend at first, and then 
increased with time. The leachate generated by MSW itself can saturated 16.7% of the landfill, so 
when designing the seepage control system and landfill leachate treatment system, this must be fully 
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considered. The soluble and degradable organic pollutants in solid phase and aqueous phase presented 
an increasing trend at first and then decreased with time, respectively, due to the release of pollutants 
from the solid phase. The peak value of the latter could reach 30 kg-m . The microorganisms in the 
solid phase and aqueous phase presented an increasing trend with time. 
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